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Abstract 

The abundant iron sulfide mineral pyrite has been shown to catalytically produce hydrogen peroxide (H 2 0 2 ) and 
hydroxyl radical ("OH) in slurries of oxygenated water. Understanding the formation and fate of these reactive 
oxygen species is important to biological and ecological systems as exposure can lead to deleterious health effects, 
but also environmental engineering during the optimization of remediation approaches for possible treatment of 
contaminated waste streams. This study presents the use of the amino acid phenylalanine (Phe) to monitor the 
kinetics of pyrite-induced 'OH formation through rates of hydroxylation forming three isomers of tyrosine (Tyr) - 
ortho-, meto-, and poro-Jyr. Results indicate that about 50% of the Phe loss results in Tyr formation, and that these 
products further react with "OH at rates comparable to Phe. The overall loss of Phe appeared to be pseudo first- 
order in [Phe] as a function of time, but for the first time it is shown that initial rates were much less than first- 
order as a function of initial substrate concentration, [Phe] D . These results can be rationalized by considering that 
the effective concentration of "OH in solution is lower at a higher level of reactant and that an increasing fraction 
of "OH is consumed by Phe-degradation products as a function of time. A simplified first-order model was created 
to describe Phe loss in pyrite slurries which incorporates the [Phe] D , a first-order dependence on pyrite surface area, 
the assumption that all Phe degradation products compete equally for the limited supply of highly reactive "OH, 
and a flux that is related to the release of H 2 0 2 from the pyrite surface (a result of the incomplete reduction of 
oxygen at the pyrite surface). An empirically derived rate constant, K pyr , was introduced to describe a variable "OH- 
reactivity for different batches of pyrite. Both the simplified first-order kinetic model, and a more detailed numerical 
simulation, yielded results that compare well to the observed kinetic data describing the effects of variations in 
concentrations of both initial Phe and pyrite. This work supports the use of Phe as a useful probe to assess the 
formation of "OH in the presence of pyrite, and its possible utility for similar applications with other minerals. 



Background 

Reactive oxygen species (ROS) are highly reactive com- 
pounds that have been studied extensively in biological 
and environmental systems and have been linked to 
numerous human health issues, including Parkinson's 
disease [1] and lung cancer [2]. In addition to forming 
naturally in cells as a function of respiration [3] and in 
the atmosphere [4], recent studies have observed that 
ROS can form at the surface of some minerals in water 
[5-7]. In particular, pyrite (FeS 2 ) has been shown to be 
efficient at forming hydrogen peroxide (H 2 0 2 ) and 
hydroxyl radical (*OH) in the presence of oxygenated 
solutions [6,8]. Of significant interest in pyrite slurries is 
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'OH, as it is transient and will rapidly react with any 
organic compound. The abundance of naturally-occur- 
ring pyrite in environments such as coal mines, where 
fine dust particles are frequently inhaled by workers, 
introduces the potential risk of human exposure to 
mineral-induced ROS. As an abundant mineral in many 
sediments and geological deposits, pyrite may also play 
a role in transformations of natural organic matter. 
Finally, the potential usefulness of pyrite as a tool for 
remediation of wastewater is also being considered as 
engineers continue to look for new methods of remov- 
ing anthropogenic compounds before discharging efflu- 
ent back into the environment. 

Several mechanisms for surface-derived ROS in sys- 
tems containing pyrite have been proposed that include 
iron-catalyzed, electron-transfer reactions. Schoonen 
et al. [6] suggests molecular oxygen is reduced to form 
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H 2 0 2 by reacting with iron-II (Fe(II)) sites at the pyrite 
surface. H 2 0 2 may either remain adsorbed on the sur- 
face and be further reduced to OH (Equation 1), or des- 
orb and undergo Fenton chemistry with dissolved 
ferrous iron (Fe 2+ ) to form OH in solution (Equation 
2); with pyrite dissolution acting as the source of Fe 2+ . 

Fe(II) (5) + H 2 0 2M -> 

Fe(III) (5) + OH" + -OH^orM 

Fe 2+ M + H 2 0 2M -> 
Fe 3+ H) + OH" + -OH H) 

The proportion of H 2 0 2 that reacts to form 'OH on 
the pyrite surface versus in solution is not known, 
although it has been proposed that the production of * 
OH and subsequent reactions with organic compounds 
occur primarily in the aqueous phase [6]. Alternatively, 
earlier work by Borda et al [8] hypothesizes that OH 
can be formed directly from the oxidation of water by 
Fe(IV) defect sites on the pyrite surface. Although more 
recent studies suggest this pathway is much less impor- 
tant in the presence of dissolved oxygen, it remains a 
potential 'OH source. 

There have been a number of studies observing inter- 
actions with OH and organic reactants in biological and 
ecological systems. Because direct measurement of 'OH 
is not possible, a number of probes have been devel- 
oped. In some studies *OH-specific phenyl-hydroxylation 
products of aromatic substrates have been employed [9]. 
In the case of pyrite slurries, probes that have been 
applied to determine 'OH-flux include electron-spin 
resonance analysis with molecular traps [10] and fluoro- 
genic probes such as 3'-(^-aminophenyl) fluorescein 
(APF) [11,12]. Cohn et al. [11] recently adapted the 
in-vitro APF method for use in pyrite slurries to quan- 
tify 'OH. However, in order to measure 'OH with APF, 
solution chemistry (e.g. pH) must be strictly controlled 
to prevent interference in fluorescence monitoring. Pyr- 
ite-mediated formation of *OH has also been implicated 
in RNA strand shortening [13] and oxidation of the 
nucleobase adenine to 8-oxoadenine [14]. Recent studies 
have shown that trichloroethylene (TCE) and its reac- 
tion products are degraded by pyrite (ultimately produ- 
cing carbon dioxide (C0 2 )), implicating pyrite-derived ' 
OH as the main oxidant [15,16]. Contributions such as 
this have led to an interest in the potential use of pyrite 
in engineered systems to facilitate remediation of 
organic chemicals in waste stream. However, the TCE 
degradation products identified are not necessarily 
specific to ROS reactions, and alternative pathways (e.g. 
microbial) may yield similar products in the environ- 
ment [17,18]. Thus there is a need for a probe to study 



reactions in unconstrained mineral systems that can be 
monitored over time with 'OH-specific products, and is 
relevant to address a range of human health and envir- 
onmental concerns. 

The fate of phenylalanine (Phe) and its degradation 
products were investigated in this work as a potential 
probe to examine pyrite-mediated 'OH reactions. As a 
naturally occurring amino acid, Phe has been shown to 
undergo *OH-specific phenyl hydroxylation reactions to 
form ortho-, meta-, and para-tyrosine (o-Tyr, m-Tyr, 
and p-Tyr) [19] as depicted in Figure 1. The o-and m- 
Tyr isomers have been used to monitor oxidative stress 
both in-vitro and in-vivo [20] as they do not form dur- 
ing normal biological processes. The Tyr isomers are 
stable enough to be measured in urine and proteins 
[21-24]. Analysis of Phe and the Tyr isomers can be 
observed simultaneously and at low levels (tens of nano- 
molar) with HPLC-MS methods, providing usefulness as 
a biologically-relevant probe. Additionally, isomers of 
Tyr can oxidize in the presence of OH to various iso- 
mers of dihydroxylphenylalanine (DOPA) (Figure 1), 
which may provide additional insight into the fate of 
Phe in mineral slurries as only 3,4-DOPA is naturally 
formed as the primary product through biological trans- 
formation of j^-Tyr [25,26]. 

This work describes the development of a sensitive 
HPLC-MS based method to evaluate the use of Phe and 
its reaction products as a way to probe 'OH-specific 
reactions involving pyrite in aqueous solution. The 
method was applied to study the kinetics of observed 
reactions as a function of both pyrite loading and con- 
centration of Phe in order to gain additional insights 
into processes that control reaction rates. As had been 
described in earlier studies with other organic com- 
pounds [6,16], the loss of Phe in an individual experi- 
ment could be adequately represented as pseudo first- 
order as a function of time. However, the effect of initial 
concentration of reactant had not been examined, and 
an unanticipated dependence on initial concentration of 
Phe ([Phe] G ) was observed. A conceptual model is pre- 
sented using a number of simplifying assumptions that 
incorporates the combined effects of pyrite loading and 
[Phe] Q on 'OH levels in bulk solution (e.g. competition 
reactions between Phe and its degradation products), 
allowing for interpretative comparisons between pre- 
dicted and observed data for both Phe and Tyr. 

Experimental Methods 

Materials 

Pyrite from Huanzala, Peru (Wards Natural Science, 
Rochester, NY) was ground and sieved to a range of 38 
- 63 [im with a surface area of 1.25 m 2 /g as per Cohn, 
et al. [14]. Phenylalanine, ^ra-tyrosine, or£/zo-tyrosine, 
3,4-dihydroxyphenylalanine, and beta-mercaptoethanol 
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(P-ME), all 99% ACS-grade or better, were obtained 
from Alfa Aesar (Ward Hill, MA). meta-Tyrosine was 
obtained from TCI America (Portland, OR). Methanol 
was GC 2 -grade from Burdick & Jackson (Morristown, 
NJ). Formic acid was ACS-grade from EMD Chemicals 
(Gibbstown, NJ). All water used for cleaning, standards, 
reactions, dilutions, and chromatography work was puri- 
fied with a Milli-Q filtration system (Millipore Corpora- 
tion, Billerica, MA) to a resistivity of 18.3 MO/cm. 

Phenylalanine degradation Experiments 

Fresh samples of crushed pyrite were treated prior to 
each experiment with a nitrogen-purged solution of 
hydrochloric acid to remove surface oxidation and 
rinsed in a glove box with nitrogen-purged water, as per 
Cohn et al. [14]. The pyrite was kept in a sealed vial 
and removed from the glove box no-more than several 
hours prior to the experiment. Aqueous stock solutions 
of Phe (or Tyr) were prepared under ambient room con- 
ditions, diluted to desired initial concentrations, and 
added to pyrite in 15 or 50 mL disposable centrifuge 
tubes to initiate the reaction (leaving 25 - 30% of the 
volume as oxic headspace). Total volume of each reac- 
tion mixture was always at least 4-times that of the 
combined volume of aliquots sequentially removed 



during kinetic studies. Tubes were then immediately set 
to rotate end-over-end on a carousel at a constant 24 
rotations-per-minute (the minimum rate found to fully 
suspend the slurry) and covered with aluminum foil to 
prevent light exposure. All experiments were conducted 
at room temperature of 25° ± 3° Celsius. 

At predetermined time points, tubes were briefly 
removed from the carrousel, briefly vortexed, and a 300 
uL aliquot was removed with a 1 mL Eppendorf® pip- 
ette at the same pyrite-to-water ratio as the sample 
(determined by a mass-balance test of repetitive sam- 
pling). Samples, including controls, were then quenched 
with 5 uL of P-ME (for a concentration of 234 mM) 
and filtered with 0.22 urn nylon Costar® centrifuge vial 
filters (Corning Life Sciences, Lowell, MA). Portions of 
the filtrate were then diluted (to different extents based 
on their initial concentrations) with water directly in a 2 
mL HPLC vial. Methanol (5%) and formic acid (0.5%) 
were also added to each vial to match initial mobile 
phase conditions. 

Each experiment was conducted independently and 
consisted of sets of incubations that ran concurrently 
for a predetermined length of time. Seven individual 
experiment sets (A - G) were conducted for this study. 
Table 1 lists the conditions for each experiment 
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Table 1 Experimental Design 
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Experiments A - G listed above with corresponding [Phe] Q , initial pyrite loading, incubation time of the experiment, and pH data. The pH was not monitored in 
every reaction vessel for each experiment as it was not a variable as all (Exp. B, C, F, and G) or multiple vials (Exp. D) had the same pyrite loadings. Reaction 
vessels in experiment G contained only one reactant each - Phe, o-Tyr, m-Tyr, or p-Tyr. Pyrite was added to the respective dilutions of Phe (or Tyr) to initiate the 
reaction. 



including incubation time, [Phe] Q levels, and pyrite load- 
ings. For the majority of the experiments conducted, 
pyrite was added at 100 g/L levels. The pH was moni- 
tored during three sets of experiments (Table 1) and 
was observed to drop rapidly within the first several 
minutes of incubation, after which it remained in a nar- 
row range over the time-course of the experiment (pH 
2.6 - 2.9) when pyrite loading was >50 g/L. When less 
pyrite was added, the final pH of these unbuffered 



solutions was higher (pH 4.2 - 5.5; Table 1), increasing 
as pyrite loadings decreased from 25 to 5 g/L. 

Instrumental Analysis 

Instrumental analyses were conducted on a Waters Cor- 
poration (Milford, MA) Alliance® 2695 HPLC coupled 
to a Waters Corporation Micromass LCT Time-of-Flight 
Mass Spectrometer (ToF-MS). A Phenomenex (Tor- 
rance, CA) Luna® C18(2) HPLC column 3 mm x 
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250 mm with 5 [im particle size was heated to 40°C and 
run with a gradient containing methanol (Solvent A) 
and 10 (iM ammonium formate/formic acid in water 
(pH 3.5) (Solvent B). Total run time was 16 minutes 
with a gradient of: 10% to 70% Solvent A over 8 min- 
utes; 70% to 10% in 4 minutes; followed by a re-equili- 
bration time of 4 minutes. A solution of Leucine 
enkephalin (Sigma Aldrich, St. Louis, MO) was injected 
post-column, generally at 1 - 3 (iL/min, for internal 
mass calibration. Mass spectrometer parameters were 
2800 V for the capillary voltage in positive-ion, electro- 
spray (ESI+) mode with cone and extraction voltages set 
at 20 V and 3 V respectively. Calibration standards con- 
taining Phe, o-, m-, j?-Tyr, and sometimes DOPA (with 
concentrations ranging from 100 nM to 5 (iM) were 
run with each sample set with method detection limits 
for Phe and Tyr of 50 nM (as determined by a signal- 
to-noise ratio of 3:1). Accurate mass measurements of 
analytes in both standards and samples were calculated 
to be within 2 mDa of the actual M+H + mass with 
spectral resolutions between 5000 - 6000 for all 
experiments. 

Controls 

Excess P-ME was added for controls in several experi- 
ments as a quenching reagent; the reported second- 



order rate constant for P-ME and *OH is 6.9 x 10 9 M" 1 
s" 1 [27], similar to 6.5 x 10 9 M" 1 s" 1 for Phe and 'OH 
[28]. No loss of Phe from solution was observed when 
0.1% (14.2 mM) P-ME was added to pyrite slurries 
(Figure 2). This provides evidence that the observed 
loss of solution-phase Phe was not due to adsorption 
of Phe to pyrite. Additional control incubations of Phe 
in water without pyrite accompanied each experiment 
and always resulted in complete recoveries, indicating 
that there were no other losses of Phe (including enzy- 
matic reactions in non-sterile media) occurring in this 
study. 

Production of H 2 0 2 has been measured in prior stu- 
dies with similarly prepared pyrite slurries [29]. No loss 
of Phe over 72 hours was observed in pyrite-free solu- 
tions of H 2 0 2 added at a wide range of concentrations 
(10 (iM, 100 (iM, and 1000 jiM) (Figure 2); nor was Tyr 
production observed. The H 2 0 2 level shown in Figure 2 
is far greater than the range found to be produced in 
160 g/L pyrite slurries, which varied from undetectable 
in the absence of the iron chelator ethylenediaminete- 
traacetate (EDTA) and up to 26 (iM in the presence of 
EDTA [29]. These observations indicate that there was 
no direct reaction between Phe and H 2 0 2 , no measur- 
able losses to adsorption of Phe to pyrite, and no micro- 
bial degradation (which would have resulted in a 




Time (hr) 

Figure 2 Degradation of various Phe concentrations in the presence of 100 g/L pyrite (Exp. A) Phe concentrations through time are fit 
well with exponential regressions. Controls of 307 uM Phe are also shown without pyrite in pure water and 1 12 uM Phe in 1000 uM H 2 0 2 ; as 
well as with pyrite quenched with 0.1% (14.2 mM) p-ME. 
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coincidental loss of Phe and preferential formation of p- 
Tyr). 

Results 

Kinetics of Phenylalanine loss 

The kinetics of Phe degradation in aqueous suspensions 
of pyrite were highly dependent on the [Phe] Q and the 
amount of pyrite mineral surface in solution. The 
observed effects of the [Phe] G proved to be interesting 
and were examined in four experiments conducted 
under different conditions (Exp. A - D; Table 1). The 
results from these studies are illustrated in Figures 2 
through 4. In each experiment the timescale for the dis- 
appearance of Phe was seen to increase with increasing 
[Phe] G . As shown in studies conducted with pyrite and 
TCE [16] or adenine [6], loss of Phe as a function of 
time is described well as pseudo first-order (Equation 3); 
however, prior studies did not examine the effect of 
initial concentrations of Phe and pyrite loadings. 



Where k* is the pseudo first-order rate constant 
derived through a fit of the data between the start of 
the reaction and when 90% of the Phe was lost (when 
reaction times permitted) (see fits in Figures 2 and 
3). Calculated k f values were inversely related to 
[Phe] Q as illustrated in Figure 2 (Exp. A) with values 
of 0.45 hr \ 0.13 hr 1 , and 0.04 hr 1 for initial con- 
centrations 26.6 uM, 103 uM, and 307 uM, respec- 
tively. Observed data was well described by Equation 
3 (Figure 1) with correlation (R 2 = 0.972 ± 0.031) 
provided in Table 2. 

Figure 2 also illustrates other important observations 
from this study. Experiment A was one of several 
experiments where a high mass loading of pyrite was 
added (100 g/L). At the highest [Phe] G , 307 uM, degra- 
dation continues for as long as 120 hours, consistent 
with evidence for long-term production of H 2 0 2 and * 
OH obtained in prior studies with pyrite conducted for 
up to several weeks under similar conditions [14]. Also, 
initial rates, R 0 (uM/hr), varied non-linearly as a func- 
tion of [Phe] Q . Values of R Q were estimated to be 4.1 
uM/hr, 9.6 uM/hr, and 10.3 uM/hr for [Phe] G of 26.6 
uM, 103 uM, and 307 uM respectively (Figure 2). R Q 
was estimated by linear regression of data collected at t 
= 0 through the time at which [Phe] is approximately 
80% of initial values. 

A similar dependence of R Q on [Phe] G was also 
observed in experiments A - D, even at different pyrite 
loadings (Exp. D) (Table 2). In experiment D, incuba- 
tions on a shorter timescale (12 hours) were monitored 
at a greater frequency for better resolution of Phe loss 
(Figure 3). The loss of Phe was determined under 



conditions where both [Phe] G (11 - 286 uM) and pyrite 
loading (25 - 100 g/L) were varied (the effect of pyrite 
loading will be addressed later in this section). 

Figure 4 shows the same relationship between the 
initial rates for each incubation from experiments A - D 
as a function of [Phe] Q when normalized to a maximum 
rate determined for each respective experiment (R Q 
/i? max ) are combined (all conducted at 100 g/L pyrite 
loading). This normalization of R 0 accounted for differ- 
ences between /? max values determined in experiments 
conducted on different days (Table 2). While pyrite 
samples were all derived from the same homogenized 
batch of ground and sieved mineral, the explanation for 
variability in the reactivity of pyrite was not determined, 
but hypothesized to have resulted from subtle differ- 
ences in acid pretreatment conditions that affected the 
abundance of Fe(III) surface sites. As further seen 
below, when results of different experiments with differ- 
ent samples of pyrite were combined, measured observa- 
tions could be described by the same mechanistically- 
based kinetic model. 

Rates increased proportionally with increasing concen- 
tration below 30 uM and clearly plateau at higher [Phe] 
0 . Estimates of R 0 at lower [Phe] G were based on few 
data points where Phe had already degraded appreciably, 
and thus estimates of R 0 are more uncertain and under- 
estimate the true initial rate to a greater extent. Still, the 
results in Figure 4 follow a hyperbolic relationship, 
which has often been interpreted using the Langmuir- 
Hinshelwood (L-H) equation when describing rates of 
reactant loss with catalysts (e.g. UV-irradiated titanium 
dioxide (Ti0 2 )) that generate 'OH at the metal surface 
[30,31]. Alternatively, it is proposed here that observed 
kinetics of the Phe data are more likely due to changes 
in ['OH](^) controlled simultaneously by the rate of 
H 2 0 2 production at the pyrite surface and OH interac- 
tion with Phe and its products (a detailed explanation 
can be found in the Discussion section). 

The effects of varying pyrite loading and surface area 
(which affect OH production) on Phe degradation rates 
are evident in Figures 3, 5, and 6. Increased pyrite load- 
ings of 25 g/L up to 100 g/L resulted in proportional 
increases in R 0 and k' at the same [Phe] G (Figure 3; 
Table 2). Figure 5 illustrates that rates (as R Q /i? max ) 
versus pyrite loading for several different [Phe] G have 
first-order dependence on pyrite loading. Although the 
effect of pH on Phe degradation was not studied for this 
work, it is noteworthy to point out that R Q determined 
at the lowest pyrite loading (pH = 5.5 at 5 g/L, Table 1 
(Exp. E)) was easily measured and still linearly corre- 
lated with those determined at higher pyrite loadings 
(pH = 2.9 and 2.6 at 50 and 100 g/L pyrite, Table 1 
(Exp. E)). Interestingly, the initial rates of H 2 0 2 forma- 
tion in similar systems have been observed to be 
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Initial Phe concentration {[Phe] 0 ): 
10 n30nM 100 nM x 300^M 



Figure 3 Concentration of Phe over time as a function of pyrite loadings (Exp. D). Data is normalized to [Phe] 0 For each pyrite loading, 
incubations of [Phe] Q from 10 uJVl to 300 uM were monitored for 12 hours. Data points for Phe are compared with fits to exponential 
regressions corresponding to a first-order model. 



relatively unaffected over a range of pHs [6]. Because 
many steps leading to the formation and fate of H 2 0 2 , * 
OH, and ferrous iron have the potential to be pH- 
dependent, more tests are needed to better understand 
the catalytic properties of pyrite as a function of pH. 

Figure 6 illustrates the relationship between Phe loss 
and the surface area of pyrite in solution (Exp. F). At 
the same pyrite loading of 10 g/L, R 0 increased from 



0.88 uM/hr to 9.4 uM/hr when the surface area was 
increased from 1.25 m 2 /g (38 - 63 urn) to 14 m 2 /g (< 38 
urn fraction) respectively (Table 2), using pyrite that 
passed through the 38- urn sieve following initial parti- 
cle-size separation. The difference in Phe loss between 
the two incubations equates to approximately a 10-fold 
increase in estimated R 0 for a corresponding 11-fold 
increase in pyrite surface area and thus is also consistent 
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Table 2 Kinetic values for experimental sets of data for this study 
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0.559 


0.9 


1.2 


1.38 


0.9915 




93.9 


25 


4.15 


0.041 


16.8 


16.9 


1.01 


0.9907 




91.9 


50 


8.65 
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6.45 
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27.7 


1.19 


0.9942 




286 
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0.047 


15.1 


14.7 


0.98 


0.9916 


E 


8.7 


5 
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0.017 
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0.8834 
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1 0.4 


50 


0.84 


0.1 30 


6.6 
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100 
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98.8 


50 


15.1 


0.079 


5.5 


8.8 
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[0-Tyr] o 94.0 


50 
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5.5 
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0.9941 




[m-Tyr] Q 97.6 


50 


21.1 


0.166 


3 


4.2 


1.39 


0.9945 




[p-Tyr] G 97.7 


50 


21.2 


0.163 


3 


4.3 


1.42 


0.9922 












Average R' 


1 (sets A, B, D, 


E, F) 


0.97 ± 0.03 


Data for R Q determined where approximately 20% of initial Phe was degraded, k' was estimated from time points up to 90% of Phe loss when possible, and 
varies between experiment sets due to variability in pyrite activity relative to OH production. Set C was a single-time point experiment and thus the data do not 
yield relevant k' values or half-lives. + = R Q was estimated as Phe loss was greater than 30% by the first time point. 



with first- order kinetics in pyrite surface area. Addition- 
ally, Figure 6 shows the Phe loss continuing for 10 days, 
twice the timescale in experiment A. 

Tyrosine production and loss resulting from 
phenylalanine degradation 

The hydroxylation products of Phe were readily mea- 
sured as o-, m-, and p-Tyv in these pyrite incubation stu- 
dies (Figure 7). The formation of the three isomers of 



Tyr has been attributed to OH-specific reactions, and 
was not observed to occur with H 2 0 2 . The production of 
the sum of the three Tyr isomers (STyr) initially occurs 
at rates that correlate with changes in initial rates of loss 
of Phe. Figure 7 illustrates the [£Tyr] through time in 
experiment A, corresponding with losses of Phe shown in 
Figure 2. Similar to R Q for Phe results, the initial increase 
in Tyr levels (from 2.0 to 2.5 (iM/hr) is relatively inde- 
pendent of [Phe] G , contrary to expectations of reactions 
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Figure 4 Initial degradation rates of Phe as a function of [Phe] Q . Shown is R Q normalized to maximum initial degradation rates (/? max ) 
each experiment set (hypothesized to vary due to subtle differences in pretreatment of pyrite). Pyrite loadings of 100 g/L were used in all 
experiments. /? max for each experiment set: A is 10.3 uM/hr, B is 4.69 uM/hr, C is 16.3 uM/hr, D is 18.9 uM/hr. Results illustrate the lack of first- 
order dependence of degradation as a function on the disappearance of Phe. 



for 



that are first-order in reactants. As the reactions proceed, 
ETyr levels peak then decrease over time. This would be 
expected as the production rate of Tyr drops due to less 
production from decreasing [Phe], combined with 
expected concurrent reactions of Phe and Tyr competing 
for the same available 'OH pool 

The highest levels of [ETyr] always peaked at time 
points occurring near to or just beyond the half-lives of 
Phe, and within a relatively narrow percentage of 
[ETyr]/[Phe] 0 , typically between 12% and 18% (Figure 
8). There is especially good agreement when [STyr] ver- 
sus Phe half-life is plotted for incubations with the same 
sample of pyrite. Results agree with predictions from the 
quantitative model developed and presented below. 

Ratios of the three Tyr isomers were observed to be the 
same among all experiments, and thus consistent with a 
common mechanism for conversion of Phe to Tyr inde- 
pendent of [Phe] G or pyrite loadings. Random hydroxyla- 
tion of the five aromatic reaction sites on Phe would result 
in a ratio of 2:2:1 for o-, m-, and ^-Tyr isomers, respec- 
tively. Illustrative of other experiments from this study, 
data normalized to [ETyr] shows that o-, m-, and p-Tyr 
were relatively constant over time in experiments with 
variable pyrite loading with ratios of 0.40 (± 0.02) : 0.28 (± 



0.02) : 0.32 (± 0.02) respectively (see Additional File 1, Fig- 
ure SI). These results show that Tyr products are not 
formed at equal-molar concentrations as hypothesized 
elsewhere [19] or at a 2:2:1 ratio. Instead, the electrophilic 
attack of the *OH appears directed to the para position 
over the ortho and meta positions. In experiments where ' 
OH was produced via pulse-radiolysis, ortho- and para- 
directed hydroxylation of toluene resulted in o-, m-, p-cre- 
sol ratios of 0.48 : 0.23 : 0.29, respectively [32]. More 
importantly, the Tyr-isomer ratios in pyrite slurries agree 
with those observed in homogeneous solutions containing 
variable levels of Fenton reagents (H 2 0 2 and Fe 2+ ) (data 
not shown). This finding is consistent with mechanisms 
where hydroxylation of Phe occurs in the aqueous phase, 
rather than a surface reaction where molecular orienta- 
tions of adsorbed species may lead to changes in reaction 
product yields. 

The degradation of individual Tyr isomers was com- 
pared to that of Phe in experiment G (Figure 9; Table 
2). Reaction rates for m- and />-Tyr were very similar to 
each other with an estimated k } of 0.17 hr 1 , which were 
faster than those determined for o-Tyr and Phe 
(approximately 0.08 hr" 1 ). Estimates for R Q for each 
reactant were closer; 22 [iM/hr for m- and p- Tyr and 
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Time (hr) 

Figure 7 Concentration of Uyr relative to [Phe] Q . The time-course of Tyr in experiment A, corresponding to the loss of Phe, illustrated in 
Figure 2. Accumulation reaches a maximum of 15% ± 3% of the initial Phe concentration at time points immediately following the half-life of 
Phe, after which degradation of Tyr exceeds its production. 



15 (iM/hr for o-Tyr and Phe. These modest differences 
in rates between Phe and Tyr are consistent with pub- 
lished second-order rate constants with 'OH of 1.3 x 



10 10 M" 



s" 1 for p-Tyr and 6.5 



10 9 M" 



s" for Phe 



(solution at pH = 2 for each) determined by pulse radi- 
olysis [28]. Thus Tyr can be expected to compete effec- 
tively with Phe with an average rate (for the 3 isomers) 
that is approximately 1.3 to 1.7 times faster. 
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Figure 8 Time-dependence of ETyr relative to [Phe] Q across all conditions of variable Phe and pyrite levels. The [ETyr] through time 
relative to the [Phe] Q peaks within a narrow range (shown in the boxed area) corresponding to similar timescales relative to Phe loss. The 
differences that are apparent in the kinetics of Tyr decay between experiments are due to variations in pyrite reactivity, but predicted by model 
calculations show below. 
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The production of DOPA as a phenyl-hydroxylation 
product of Tyr was also observed in this study. DOPA 
was identified by correspondence of HPLC retention 
with an authentic standard of one of the possible iso- 
mers (3,4-DOPA) which eluted before p-Tyr, and 
further confirmed by accurate mass measurements by 
Tof-MS that were within 2 mDa of the actual mass. 
This finding is consistent with production of DOPA for 
reactions of Tyr and 'OH in prior work [26]. No other 
peaks corresponding to the mass of DOPA were identi- 
fied. Due to the very low concentrations of DOPA 
(when detected) in Phe experiments at the dilution- 
levels injected, it was not routinely monitored. It is 
noteworthy that in an experiment with p-Tyr, DOPA 
was only measured at a maximum concentration around 
2% relative to [p-Tyr] Q . Although a comprehensive 
HPLC-MS method to chromatographically separate, 
indentify, and quantify the 6 possible isomers of DOPA 
was not conducted, it appears that Tyr-to-DOPA may 
not be as sensitive of a probe as Phe-to-Tyr in monitor- 
ing *OH-specific reactions mediated by pyrite. 

Discussion 

Phenylalanine as a hydroxyl radical probe 

A sensitive method was developed for the quantification 
of Phe and Tyr at >50 nM levels by HPLC-MS methods 



with direct aqueous injection without the need for pre- 
concentration or derivatization steps that could lead to 
more analytical uncertainty. The role of 'OH in the 
degradation of Phe is confirmed by the production of a 
characteristic composition of three isomers of Tyr - 
ortho, meta, or para. The presence of m- and o-Tyr has 
been used previously to assess the importance of OH in 
oxidative stress [20]; whereas observation of primarily p- 
Tyr is normal in biological systems. High yields, charac- 
teristic ratios, and persistence of readily-measured Tyr 
products are other traits that make Phe a good probe to 
monitor *OH production in pyrite and other mineral 
slurries. Yields of Phe to STyr conversion were esti- 
mated in this study to be between 30% and 60% by 
extrapolating initial changes in the [STyr] t /([Phe] Q - 
[Phe] t ) back to time-zero (data not shown). The stability 
in pyrite slurries of Tyr is similar to that of Phe, allow- 
ing it to be monitored at appreciable levels throughout 
the period of Phe decay. 

There are of course other molecular probes that have 
been developed that rely upon measurement of phenyl 
hydroxylation products to determine concentration 
levels of 'OH in aqueous solutions. Two of the more 
sensitive ones are production of hydroxybenzoic acids 
from benzoic acid [9] and hydroxyterephthalic acid from 
terephthalic acid [33]. In each case, the determination of 
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products by fluorescence or HPLC-fluorescence can be 
even more sensitive and require less expensive instru- 
mentation than the HPLC-MS based detection of Phe 
and Tyr described here. However, neither benzoic acid 
nor terephthalic acid are nearly as fluorescent as their 
products. Hence, in experiments with those probes only 
the hydroxylated products are typically measured, while 
the substrate - often added in excess - is routinely not 
analyzed for. When detection of these substrates is 
required, it is usually conducted with far less sensitive 
UV-based methods [34]. 

There are other potential advantages of developing 
Phe as a probe to study *OH-mediated reactions. The 
presence of Phe in cells and body fluids allows for 
meaningful interpretation of the ratios of o-, m-, p-Tyr 
in vivo [21,23]. Similarly, in contrast to other aromatic 
probes used to date, Phe is naturally present at low 
levels in all natural waters and biologically active geolo- 
gical matrices, including environments where pyrite is 
present. With methods available for determining ultra- 
trace environmentally relevant concentrations of Phe 
and Tyr isomers [35,36], it may be possible to conduct 
field studies that could shed insight into the conditions 
where 'OH-mediated reactions are occurring in the 
environment. 

The mass spectrometric methods utilized also provide 
certain advantages. With the full-spectral sensitivity of LC- 
ToF-MS, direct injection of aqueous samples offers the 
potential to identify other non-targeted reaction products 
that may have otherwise been lost in isolation or purifica- 
tion steps. For example, in this work, trace levels of DOPA 
were identified in some experiments with confirmation 
using accurate-mass measurement. Finally, in more com- 
plex matrices characteristic of biological fluids or organic 
matter-rich natural waters, the increased specificity of ana- 
lysis by ToF-MS allows for greater discrimination from 
matrix that can potentially complicate the use of some UV 
or fluorescence-based detection methods. 

Pyrite-mediated hydroxyl radical formation 

The mechanisms of OH production in pyrite slurries 
are not certain, yet appear important for understanding 
the reaction it undergoes with Phe and other organic 
compounds. For example, if *OH (fl(is ) is the primary 
source of OH involved in Phe reactions, and Phe reacts 
at the pyrite surface, then the dependence of R 0 on 
[Phe] G might be appropriately described by a L-H sur- 
face catalysis model. However, if the primary source of * 
OH is derived in solution by Fenton-like reactions, then 
there needs to be a reaction mechanism to describe the 
nonlinear response of R 0 with varying [Phe] G (Figure 4). 

Any surface catalyzed reactions would most likely 
require the production of adsorbed 'OH which may be 



formed through cathodic reduction of H 2 0 2 at the pyr- 
ite surface (Equation 1) [6], or through Fe(II)/Fe(IV) 
electron-transfer reactions with water at surface defect- 
sites [8]. However, it has been argued recently that sur- 
face-defect sites in oxygenated slurries are less likely to 
form "OH than the Fenton reaction in solution [6]. 
Reaction of Phe directly at the pyrite surface with OH 
{ads) seems unlikely to be a major contributor to its 
overall loss as Phe, Tyr, and soluble probes used in 
prior studies of pyrite mediated reactions [14] do not 
measurably adsorb to the surface (Figure 2). However, 
desorption of *OH ( ^ 5) may occur and resulting reactions 
(i.e. with Phe) could be limited to diffuse boundary 
layers when sufficient substrates are within a limited dis- 
tance from the surface [30]. Turchi and Ollis [30] 
showed that such boundary layer reactions could explain 
kinetic behavior consistent with the L-H model, provid- 
ing one possible mechanism for near surface reactions 
that cannot be completely ruled out. 

The possibility that higher Phe levels could affect the 
catalytic properties of pyrite seems unlikely as in the 
absence of observed adsorption of Phe, its surface cover- 
age would be very low, and decrease with increasing 
pyrite loading. Furthermore, there was no observed 
effect to reaction rates with the addition of excess Fe 2+ 
(see below), consistent with Phe not affecting the avail- 
ability of iron in solution needed for the Fenton 
reaction. 

Support for *OH formation via the Fenton reaction in 
the aqueous phase (Equation 2) is much stronger for a 
number of reasons. For example, dissolved H 2 0 2 can be 
observed accumulating in pyrite slurries at concentra- 
tions relevant to this and other studies (micromolar- 
range) when the Fe 2+ in solution is chelated by EDTA, 
inhibiting the Fenton reaction [6,16,29]. Ferrous iron in 
solution also acts as a catalyst in the Fenton reaction 
and has been measured elsewhere to be in excess of 
H 2 0 2 [29]. To confirm that ferrous iron is present in 
excess, experiments with 50 g/L pyrite slurries contain- 
ing excess ammonium ferrous sulfate (500 (iM) were 
incubated with 100 \iM Phe. Results indicated no differ- 
ence in initial or overall rates of Phe degradation when 
compared to a pyrite solution without iron addition (see 
Additional File 1, Figure S2). Therefore, there is suffi- 
cient Fe 2+ in pyrite slurries to suggest the OH produc- 
tion is limited to the rate of production and release of 
H 2 0 2 , and thus by extension, the production rate of 
H 2 0 2 at the pyrite surface affects the kinetics of Phe 
loss (and its degradation products). 

Kinetics of Phenylalanine loss 

The reaction of OH and Phe is described by the follow- 
ing second-order rate equation: 
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d[Phe] 
dt~~ 



fep he [Phe] [-OH] 



(4) 



Where A- phe is the second-order rate constant and is 
equivalent to 6.5 x 10 9 M" 1 s" 1 (as mentioned in the 
Results section). However, it is shown (e.g. in Figures 2, 3, 
4) that -d[Phe]/dt does not vary proportionally with [Phe] 
o, and plateaus at higher concentrations. We suggest that 
the most likely mechanism for these observations is that 
as [Phe] G increases, there is a near-proportional decrease 
in [*OH] in the aqueous phase. This situation contrasts 
with more typical experiments monitoring the production 
of 'OH or its reaction products, where either excess reac- 
tant is used and product determined, or in cases where the 
reactant is added at low levels in comparison to other spe- 
cies in solution that are sinks for *OH. This balance is con- 
strained by the limited production rate (or flux) of OH at 
any given time. To support this hypothesis we present a 
conceptual model that can explain most of the results 
from this study using three primary assumptions. 

Assumption 1 asserts that the OH flux is constant 
over the time course of the experiment. This is supported 
by evidence that H 2 0 2 production has been shown to 
be proportional to pyrite loading [29], and does not 
react with substrates such as adenine [14] and Phe. If 
the rate of OH formation is assumed to be constant 
and proportional to pyrite loading, it can be assigned a 
reaction coefficient that represents the pyrite reactivity 
(Equation 5), incorporating all factors related to forma- 
tion of *OH for a particular sample of pyrite. 



R OH formation = Kpyr[pvr] 



(5) 



Where K pyr has units of mol g" 1 hr" 1 (and also corre- 
sponds to the maximum rate at which Phe can be 
degraded). K pyr is determined for each pyrite sample to 
account for aforementioned variability in reaction rates 
(i.e. rate of *OH-formation) that was observed between 
different experiment sets and hypothesized to be related 
to conditioning of the pyrite prior to incubation. 

Assumption 2 states that OH is in approximate steady- 
state with a loss dominated by reaction with Phe (and its 
degradation products) and the rate of formation related to 
the loading and surface area of pyrite. The rate of decay 
for each reactant product through time is related to its 
individual rate constant and concentration, and with 'OH 
formation assumed constant in Equation 5, steady-state 
requires the loss of OH be constant. Thus the combined 
loss of total reactants must also remain constant. 



d[OH] 



Kpyrlpyr]- 



dt pyrl 
•OH](fe Phe [Phe] + fe T yr[Tyr] + 

E(fei[i])) = 0 



(6) 



Where Ar Tyr and k x are the second-order rate constants 
for the reaction of OH with Tyr and other degradation 
intermediates, collectively S U respectively. 

Equation 6 represents a balance of production and loss 
of 'OH, but also underscores that many degradation pro- 
ducts of Phe that are likely important for describing the 
kinetics are unaccounted for. The degradation intermedi- 
ates compete with Phe, Tyr, and each other for reaction 
with the limited OH until completely oxidized to carbon 
dioxide (C0 2 ). Therefore, the third simplifying assumption 
states that the sum of all oxidizable reactants is equal to 
[Phe] 0 (i.e. [Phe] t + [Tyr] t + Y[i] t « [Phe] G ) for most of the 
time course of the experiments (i.e. where Phe was still 
readily detectable (>10% Phe remaining)). In addition, it is 
also assumed that the second-order rate constants of those 
products with 'OH are all relatively high and similar to 
each other (i.e. k Phe « k Tyr « ki). Similar approaches have 
been used to simplify complex competitive rate equations 
in other studies involving surface-catalyzed *OH-mediated 
reactions [31]. Therefore, a simplified approximation of 
Equation 6 is represented in Equation 7. 



d[OH] 



%r[pyr]- 



dt ^ pyn 
fe Phe [Phe]o[-OH] = 0 

Which can also be rewritten in terms of [*OH] 



•OH] 



K pyr[pvr] 
fc Phe [Phe] 0 



(7) 



(8) 



With substitution of Equation 8 into Equation 4, Phe 
loss can now be described with a rate expression that is 
a function of both pyrite loading and [Phe] G . 



d[Phe] K pyr [pyr] 
dt " [Phe] G 



[Phe] 



(9) 



Note that in Equation 9, K pyr [pyr]/[Phe] G is analogous 
to the experimentally determined pseudo first-order rate 
constant k } defined in Equation 3, where k y was found 
to be proportional to pyrite loading and inversely related 
to [Phe] G (at [Phe] > 30uM). 



K pyr[pvr] 
[Phe] G 



(10) 



Therefore, K pyr can be derived experimentally for a 
given pyrite sample and allows for the modeling of Phe 
and Tyr concentrations through time. 

Modeling Phenylalanine and Tyrosine concentrations 
through time 



Determination of K. 



pyr 



To account for the differences in reactivity among pyrite 
samples, a unique value for K pyr was determined for 
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experiments A, B, D, and E using all kinetic data col- 
lected. As described in the Results section, the temporal 
loss of Phe generally fits an exponential function well 
where k y is the slope of In [Phe] versus t. As per Equa- 
tion 10, K pyr is then derived from the slope of k' versus 
[pyr]/[Phe] Q . The measured data were described by lin- 
ear correlations between k' and [pyr]/[Phe] Q (see Addi- 
tional File 1, Figure S3), and an averaged K pyr was 
estimated for each experiment (A, B, D, and E) by set- 
ting the y-intercepts to zero. 
Phenylalanine loss 

Integration of Equation 9 allows for the modeling of Phe 
loss through time. 



/ Kp y r[pyr] \ 



[Phe] t = [Phe] 0 



[Phe] Q 



V 



(11) 



/ 



Observed Phe data from Experiments A, B, D and E 
along with predicted values (modeled with Equation 11) 
are plotted in Figure 10. Comparisons of predicted 



versus observed [Phe] t was fit by a linear regression 
model with a slope of 0.985 with an R 2 of 0.996. 
Tyrosine production and loss 

The kinetics of Tyr formation and loss are governed by 
consecutive reactions, and using the assumptions 
described above, the rate of change can be represented 
in Equation 12. 



dJTyr] 
dt 



Kpyr[pyr]- 

afcphe[Phe] - fer yr [Tyr] 
fcphe[Phe] + fc Ty r[Tyr] + £(fci[i])) 



(12) 



Where a is the fraction of the Phe reaction products 
corresponding to Tyr. For modeling purposes, a was 
estimated to be 0.5, which was within the range of 
yields determined from initial rate data. Again, 
assumptions made in this model consolidate the com- 
petitive factor that affect Phe and other products that 
react with OH in the slurry. Thus Equation 12 can be 
simplified to: 
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Figure 10 Modeled versus observed Phe concentrations for 4 different sets of experiments. Calculated concentrations of Phe were 
derived from Equation 10 using the K pyr calculated for experiment sets A, B, D, and E. The model values fit the observed data well despite the 
observed variability in K pyr among batches of pyrite. 
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d[Tyr] K pyr [pyr] 



(a [Phe] - [Tyr]) 



dt [Phe] D 
Integration of Equation 13 yields: 

/ fepyr[pyr] 

[Tyr], = a K pyr [pyr]t e l phe J„ 



(13) 



(14) 



Equation 14, based on data for Phe loss, provided pre- 
dictions of [Tyr] t that agreed with measured values from 
Experiments A, B, D, and E reasonably well (Figure 11). 
A linear regression fit of points in Figure 11 has a slope 
of 1.18 with an R 2 of 0.888. 

This simplified model is able to predict much of the 
change in Phe and Tyr through time when the specific 
K pyr for a pyrite sample is determined. Importantly, it 
allows for most of the observed dependence of [Phe] G 
and pyrite loading to be accounted for while describing 
the apparent first-order dependence of Phe loss through 
time. Based on this analysis, there is no need to consider 
potential surface reactions between OH and Phe to 



model the hyperbolic function of R Q at higher [Phe] G 
(Figure 4), normally ascribed to surface catalytic reac- 
tions. However, it should be pointed out that this simple 
model does not describe the decrease in R 0 observed at 
lower [Phe] G as well. Although this may be related to 
the paucity of initial measurements when Phe degrades 
rapidly at lower [Phe] G , it does seem likely that initial 
rates are lower when [Phe] G is low (<30 uM). We 
hypothesize that at low [Phe] Q , there are other reactions 
that compete for 'OH with Phe and other products. 
Whether or not such reactions occur at the pyrite sur- 
face or with trace inorganic reactants in solution is 
uncertain and beyond the scope of this work. 

A more detailed numerical simulation was also devel- 
oped to compare to results from the simplified model 
presented above (Equations 11 and 14), and to provide 
greater flexibility for testing the effect of assumptions 
concerning reaction pathways and relative rates of mul- 
tiple products with OH (see Additional File 1, Numeri- 
cal simulation section). Values of second-order rate 
constants for different reaction products (i.e. Tyr and 
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Figure 11 Modeled versus observed Uyr concentrations for 4 different sets of experiments. Corresponding to the data for Phe shown in 
Figure 10, concentrations of ZTyr were derived using Equation 14, the K pyr calculated for each experimental set, and a value of 0.5 for a. Note 
that for points where values are predicted but not measured for Tyr, there was an apparent lag in Phe degradation at the first time point (see 
inset). Similarly, several instances when no Tyr is predicted to be present, low-levels were observed corresponding to later time points in the 
incubation where the vast majority of Tyr had been degraded. 
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DOPA) and the fraction of Phe-to-STyr conversions 
(reaction yield, or a) were adjusted to test basic assump- 
tions made above. The assumption of constant *OH-flux 
was not changed. Figures in the Additional File 1 show 
that calculations made with Ap he = /c Tyr = /cdopa and a 
= 0.5 (50% conversion of Phe-to-ETyr) offer the closest 
match to observed and modeled data for most 
incubations. 

Conclusion 

Loss of Phe varied in pyrite slurries at rates that were 
first-order in pyrite loading and pseudo first-order 
dependent on Phe as a function of time; whereas the 
rate of Phe loss was much less-than first-order in [Phe] 
G . The data for the loss of Phe as well as the production 
and loss of Tyr products could be described well by a 
mechanistically-based kinetic model that reconciles the 
observations concerning the initial concentration depen- 
dence of Phe. The competitive effects of degradation 
products on reactions of molecular probe have been 
included in kinetic descriptions; as such products 
become important when substrate is not added in great 
excess. Not limited to experiments with pyrite, there 
have been few studies that have considered that 
observed exponential decay of reactants may be due to 
increased competition for available *OH [31], rather 
than due to true first-order behavior. In systems where 
the flux of OH is likely the rate limiting step, observed 
kinetics will depend on whether the relative amounts of 
competing reactants change as a function of time. The 
[STyr]/[Phe] 0 and [Tyr-isomer]/[STyr] ratios were con- 
sistent throughout experiments, with total-Tyr forma- 
tion estimated to be about 50% of the Phe conversion 
via OH. The use of Phe and its *OH-specific products is 
argued to be a useful probe that should be further 
developed for the study of the mechanisms of pyrite and 
other mineral-mediated reactions, and has the potential 
to be a valuable tool for the study of 'OH reactions in a 
range of other systems with more complex matrices. 
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